Aspects of Hofava-Lifshitz cosmology 
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We review some general aspects of Ifofava-Lifshitz cosmology. Formulating it in its basic version, 
we extract the cosmological equations and we use observational data in order to constrain the pa- 
rameters of the theory. Through a phase-space analysis we extract the late-time stable solutions, and 
we show that eternal expansion, and bouncing and cyclic behavior can arise naturally. Concerning 
the effective dark energy sector we show that it can describe the phantom phase without the use of 
a phantom field. However, performing a detailed perturbation analysis, we see that Hofava-Lifshitz 
gravity in its basic version suffers from instabilities. Therefore, suitable generalizations are required 
in order for this novel theory to be a candidate for the description of nature. 
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I. INTRODUCTION 



Almost one year ago Hofava proposed a power- 
counting renormalizable theory with consistent ultra- 
violet (UV) behavior \lttM- Although presenting an in- 
frared (IR) fixed point, namely General Relativity, in the 
UV the theory exhibits an anisotropic, Lifshitz scaling 
between time and space. Due to these novel features, 
there has been a large amount of effort in examining and 
extending the properties of the theory itself [3-[Tl|. Ad- 
ditionally, application of Hofava-Lifshitz gravity as a cos- 
mological framework gives rise to Hofava-Lifshitz cosmol- 
ogy, which proves to lead to interesting behavior [l^] ■ In 
particular, one can examine specific solution subclasses 
[Tsl - fTB l. the phase-space behavior [l6l-fl8!|. the gravita- 
tional wave production (lOi^he perturbation spectrum 
[2(]| - |22i . the matter bounce [23l - l26| , the black hole proper- 
ties |27H29| , the dark energy phenomenology fSOUSSl] , the 
observational constraints on the parameters of the the- 
ory [s^ - lssj , the astrophysical phenomenology [37] , the 
thermodynamic properties jssl . Isoj etc. However, despite 
this extended research, there are still many ambiguities 
if Hofava-Lifshitz gravity is reliable and capable of a suc- 
cessful description of the gravitational background of our 
world, as well as of the cosmological behavior of the uni- 
verse liMiilHIia. 



In the present work we review the basic aspects of 
Hofava-Lifshitz cosmology. The manuscript is organized 
as follows: In section |ll] we present the simple version of 
Hofava-Lifshitz cosmology, in both its detailed-balance 
and beyond-detailed-balance version. In section IIIII we 
use observational data in order to constrain the parame- 
ters of the theory. In section ITVl we present the results of 
the phase-space analysis, in |V] we present the bouncing 
and cyclic solutions, and in IVII we extend the theory in 
order to present a more realistic dark energy phenomenol- 
ogy. In section IVlIl through a perturbation analysis, we 



discuss the instabilities of the simple versions of the the- 
ory, and thus in section IVIIII we present a healthy ex- 
tension of Hofava-Lifshitz gravity. Finally, section IIXI is 
devoted to the summary of our results. 



II. HORAVA-LIFSHITZ COSMOLOGY 

In this section we briefly review the scenario where the 
cosmological evolution is governed by the simple version 
of Hofava-Lifshitz gravity [l^] . The dynamical variables 
are the lapse and shift functions, N and Ni respectively, 
and the spatial metric gij (roman letters indicate spatial 
indices) . In terms of these fields the full metric is written 



ds^ = -N'^dt^ + g,j{dx' + N'dt){dx^ + N^dt), (1) 

and the scaling transformation of the coordinates reads 
as < -> ^^t and x'' Ix"^ . 



A. Detailed Balance 

The gravitational action is decomposed into a kinetic 
and a potential part as S'g = / dtd^Xy/gN{CK+Cv)- The 
assumption of detailed balance [3] reduces the possible 
terms in the Lagrangian, and it allows for a quantum 
inheritance principle [l|, since the {D + l)-dimensional 
theory acquires the renormalization properties of the D- 
dimensional one. Under the detailed balance condition 
the full action of Hofava-Lifshitz gravity is given by 



+ 



dtd^x^N { —{KijK'^ - XK^) 



8(3A- 1) 



- 4A 



Rii^ jRk 



-W + AR- 



(2) 
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where Kij — {glj — WiNj — WjNi) /2N is the extrinsic 
curvature and C'^ = e*J'=Vfc(i?^ - R5l/4)/^ the Cot- 
ton tensor, and the covariant derivatives are defined with 



2 



respect to the spatial metric gij. e^^^ is the totally an- 
tisymmetric unit tensor, A is a dimensionless constant 
and the variables k, w and /i are constants. Finally, we 
mention that in action ^ we have already performed 
the usual analytic continuation of the parameters and 
w of the original version of Hof ava-Lifshitz gravity, since 
such a proccdmc is required in order to obtain a realistic 
cosmology 13, 14, 28, 38] (although it could fatally affect 
the gravitational theory itself). 

In order to add the matter component we follow the hy- 
drodynamical approach of adding a cosmological stress- 
energy tensor to the gravitational field equations, by de- 
manding to recover the usual general relativity formula- 
tion in the low-energy limit jB, 16, 32]. Thus, this matter- 
tensor is a hydrodynamical approximation with pm and 
Pra (or Pm and Wm) as parameters. Similarly, one can 
additionally include the standard-model-radiation com- 
ponent, with the additional parameters pr and Wr- 

In order to investigate cosmological frameworks, we 
impose the projectability condition [5] and we use an 
FRW metric 



with 



(3) 



(4) 



where K <,—,> corresponding to open, flat, and 
closed universe respectively. By varying N and gij, we 
extract the Friedmann equations: 



H 



6(3A- 1) 
6(3A- 1) 
8(3A- l)2a2 



2 4(3A-1) 



8(3A- l)a4 ^ 8(3A- 1) 



(5) 



WmPrn + WrP' 



4(3A- 1 

16(3A- l)2a2 



8(3A-l)a4 8(3A-1) 



(6) 



where H = ^ \s the Hubble parameter. As usual, pm 
follows the standard evolution equation /9,„ -I- 3H{p,n + 
Pm) = 0, while Pr follows Pr + 3H(pr + Pr) = 0. Finally, 
concerning the dark-energy sector we can define 



PDE 



Pde 



8(3A- l)a4 8(3A~ 1) 



k2^2^2 



3k2^2a2 



8(3A-l)a4 8(3A-1)' 



(7) 



(8) 



The term proportional to a"** is the usual "dark radi- 
ation term", present in Hof ava-Lifshitz cosmology [l^ . 
while the constant term is just the explicit cosmologi- 
cal constant. Therefore, in expressions ([T]),® we have 
defined the energy density and pressure for the effective 
dark energy, which incorporates the aforementioned con- 
tributions. Note that using Q,© it is straightforward to 
show that these dark energy quantities satisfy the stan- 
dard evolution equation: pde + 'iH{pDE + Pde) — 0. 

If we require expressions ([5|) to coincide with the stan- 
dard Friedmann equations, in units where c = 1 we set 



K 



8(3A- 1)2 



1, 



(9) 



where Gcosmo is the "cosmological" Newton's constant, 
that is the one that is read from the Friedmann equa- 
tions. We mention that in theories with Lorentz invari- 
ance breaking Gcosmo does not coincide with the "grav- 



itational" Newton's constant G„ 



that is the one that 



is read from the action, unless Lorentz invariance is re- 
stored ll^. For completeness we mention that in our 
case Ggrav = K^/(327r), as it can be straightforwardly 
read from the action ([2]). Thus, it becomes obvious that 
in the IR (A 1), where Lorentz invariance is restored, 
Gcosmo and Ggrav coiucide. 



B. Beyond Detailed Balance 

The aforementioned formulation of Hofava-Lifshitz 
cosmology has been performed under the imposition of 
the detailed-balance condition. However, in the litera- 
ture there is a discussion whether this condition leads to 
reliable results or if it is able to reveal the full informa- 
tion of Hofava-Lifshitz gravity [l^ . Therefore, one needs 
to investigate also the Friedman equations in the case 
where detailed balance is relaxed. In such a case one can 
in general write [1-13, d, [13 : 



2ao 



(3A-1) 
2 

(3A^ 

0-2 



6 
K 



Pr 

6a4 



<J4K 



6a6 



3(3A- 1) a2 



(10) 



H 



(3A-1) 
3 

(3A^ 

CT2 



6 



lUrPr 



18a4 



6a6 



K 



6(3A- 1) a? 



(11) 



3 



where (Tq = and the constants ai are arbitrary 

(with (72 being negative and a4 positive). Furthermore, 
the dark-energy quantities are generahzed to 

Pde\ 



11 
6 



11 

' 6 



6a4 



6a6 



18a4 



6a6 



(12) 
(13) 



Again, it is easy to show that 

P-D-B|„o„-db + 3-ff(/3_D_ELo„_di 

Finally, if we force p^ . (|lip to coincide with the standard 
Friedmann equations, we obtain: 



+ Pdb|„)=0. (14) 



6cro 



' 87r(3A-l) 
-3(3A- 1), 



(72 = -3(3A- ij, (15) 

while in this case the "gravitational" Newton's constant 
Ggrav writes as Gg^av = 6ao/(167r). Similarly to the de- 
tailed balance case, in the IR (A = 1) Gcosmo and Ggrav 
coincide. 



III. OBSERVATIONAL CONSTRAINTS 

Having presented the cosmological equations of a uni- 
verse governed by Hofava-Lifshitz gravity, both with and 
without the detailed-balance condition, we now proceed 
to study the observational constraints on the model pa- 
rameters 1341 [35I. 



A. Constraints on Detailed-Balance scenario 

We work in the usual units suitable for observational 
comparisons, namely setting 87rGgiav = 1 (we have al- 
ready set c = 1 in order to obtain (jO])). This allows us 
to reduce the parameter space, since in this case Ggiav 

lead to: Gcosmo 



gives K — 4 and thus 

,2 A _ (3A-r" 



an ii^K = 



47r(3A-l) 

In order to proceed to the elabora- 



tion of observational data, we consider as usual the mat- 
ter (dark plus baryonic) component to be dust, that is 
Wm ~ 0, and similarly for the standard-model radiation 
we consider Wr = 1/3, where both assumptions are valid 
in the epochs in which observations focus. Therefore, 
the corresponding evolution equations give pm = Pmo/a^ 
and pr = Pro/ci^ respectively. Additionally, instead of 
the scale factor it proves convenient to use the red- 
shift z as the independent variable, which is given by 
\ + z = ao/a = 1/ a. Finally, we introduce the usual den- 
sity parameters [VLm = Pm/{3H'^), fix = ~K/{H^a'^), 
ilr = pr/{3H'^))- Inserting these relations into Fried- 
mann equation ([5]) we acquire: 

2 



KO 



(3A-1) 
(1 + -)^ 



+ %l(l + .)^]}, (16) 



where we have also introduced the dimensionless parame- 
ter uj = t5-4t , and where a 0-subscript denotes the present 

value of the corresponding quantity. Applying this rela- 
tion at present we get: 



(3A-1) 



nr. 



KO 



02 



1. 



(17) 



We remind that the term JI^q/ i^^) ^^'^ coefficient of 
the dark radiation term, which is a characteristic feature 
of the Hofava-Lifshitz gravitational background. Since 
this dark radiation component has been present also dur- 
ing the time of nucleosynthesis, it is subject to bounds 
from Big Bang Nucleosynthesis (BBN). As discussed in 
more details in the Appendix of [34] , if the upper limit on 
the total amount of dark radiation allowed during BBN 
is expressed through the parameter ANi, of the effective 
neutrino species ^43j], then we obtain the following con- 
straint: 



''''KO 



0.135A7V^r2rt. 



(18) 



In summary, the scenario at hand involves four pa- 
rameters (we fix Hq by its 7-year WMAP best-fit val- 
ues, given in Table 1 of 144]), namely VLmo, ^ko, ^ and 
AiV^, subject to constraint equations (|17p and (fT8)) . We 
marginalize over the cosmological parameters rimOj ^^fcOi 
firo and Hq. Of the four remaining parameters, only two 
are independent, and we choose A and l^N^, as our free 
parameters. Once these are chosen, and for a given choice 
of curvature, O/fo and lo are immediately fixed from the 
constraint equations. In particular, w can be determined 
by eliminating VIko from relations p7)) and (IT8| : 



w-2 sgn {nxo) VO^TsSAiViTO^ 

rimO + n^o 



3A- 1 



1 = 0. 



0.135AiV^Ort 

(19) 



VLko can then be found from lj using (|18p . 

In Fig. [l]we use a combination of observational data 
from SNIa, BAO and CMB to construct likelihood con- 
tours for the parameters rimo and AiV^ for positive cur- 
vature. Additionally, in Fig. [5] we display the likelihood- 
contours for the free parameters A vs AiV^ for posi- 
tive curvature, where all other parameters have been 
marginalized over. Finally, in Table U we summarize the 
1(7 limits on the parameter values for the detailed-balance 
scenario. 



K 




[l/Hl) A 


(87rGgrav-ffo) M 


A 




>0 


4 


(0, 1.46) 


(1.37, cx)) 


(0.98, 1.01) 


(0, 0.32) 


<0 


4 


(0, 1.46) 


(1.8,00) 


(0.97, 1.01) 


(0, 0.68) 



TABLE I: 1(7 limits on the parameter values for the detailed- 
balance scenario, for positive and negative curvature. The 
cosmological parameters f2mO, ^bo, firo and Ho have been 
marginalized over. 



4 



0.31 



, 0.29 



0.27 



0.25, 




where we have introduced the dimensionless parameters 
wi = = """l^^'"^ and UJ4 = _£4^. Addition- 

ally, we consider the combination lu^ to be positive, in 
order to ensure that the Hubble parameter is real for all 
redshifts. 

In summary, the present scenario involves the following 
parameters: the cosmological parameters Hq, flmo, flxo, 
f2f,0i f^rOj and the model parameters A, wi, W3 and UJ4. 
Similarly to the detailed-balance section these are subject 
to two constraints. The first one arises from the Friedman 
equation at 2; = 0, which leads to 



FIG. 1: (Color Online) la and 2a contour plots of Q,mO vs 
ANi, for positive curvature (K > Q), under SNIa, BAO and 
CMB observattonal data. The whtte diamond marks the best- 
fit point. The model parameters co = A/(2Hq) and Q,ko = 
—K/{Hq) are related to flmo and ANi, through fl8\) and fl9\) . 



(3A-1) 



nKo = i- (20) 



This constraint eliminates the parameter wi. The second 
one arises from BBN considerations, since, as we show in 



the Appendix of [3J|, at the time of BBN (z 
acquire [43j : 



zbbn) we 



1.06 



1.02 




ws + a;4 (1 + zbbn)^ = t^sn 



Q.135ANMro, (21) 



0.98 



0.94 



FIG. 2: (Color online) la and 2a contour plots of X vs AN^ 
for positive (K > Q) curvature, in the detailed-balance sce- 
nario, under SNIa, BAO and CMB observational data. The 
remaining parameters have been marginalized over. 



In conclusion, we see that the Hofava-Lifshitz cosmo- 
logical scenario under the detailed balance condition is 
not ruled out by observations. However, there are tight 
constraints on the model parameters. Furthermore, the 
data constrain A to roughly A = iIq 02 1'^ level, 

that is to a very narrow window around its IR value, 
while its best fit value is very close to I {Xb.f = 0.006). 



B. Constraints on Beyond-Detailed-Balance 
scenario 



In imits where STrGgi-av — 1, Ggrav gives gq — 1/3. 
Following the procedure of the previous subsection, the 
Friedmann equation (|10p can be written as 



(1 + + J^^—^ [^1 +L03{1+ Zf + (1 + ^)'] }, 



where wsmax denotes the upper limit on W3. In the 
following, we use expression (j2ip to eliminate 0-14. For 
conveniencejinstead of ^3 we define the new parameter 



We use relation (1211) to eliminate 0^4 in favor of a and 
ANi,, and treat A, a, Qko and AN^, as our free param- 
eters, marginalizing over Hq, flmo, ^bo and flro. Using 
the combined SNIa-t-CMB-t-BAO data, we construct like- 
lihood contours for different combinations of the above 
parameters. Figure [3] depicts the Icr and 2a W3 — ir^Kol 
contours, for ATV^ = 2, for positive curvature, while Fig. 
m depicts the A-variation. The approximate la limits 




iog^„(w3) 



FIG. 3: (Color Online) la and 2a contour plots 0/ logj^g(u;3) 
vs logio \Qko\ for K > and AN^ = 2.0, using SNIa, BAO 
and CMB data. The white diamond marks the best-fit point. 



on the model parameters ai are presented in Tabic [TTl 
Additionally, in Table HH] we focus on the la limits of a 
and A. 
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cro/(87rG) 


AN. 




(S-KG/Hi) ai 


(72 


(S-KGHi) (73 


cr4/(87rG) 


1/3 


0.1 


(0, 0.01) 


(4.29, 4.33) 


-6 


(0, 0.03) 


(-9.08 X 10"^^ 0) 


1/3 


0.1 


(-0.01, 0) 


(4.40, 4.45) 


-6 


(0, 0.81) 


(0, 5.66 X 10"^^) 


1/3 


2.0 


(0, 0.04) 


(4.13, 4.45) 


-6 


(0, 0.01) 


(-1.77 X 10"^", -2.62 X 10-^^) 


1/3 


2.0 


(-0.01, 0) 


(4.40, 4.45) 


-6 


(0, 0.23) 


(-2.61 X 10~^'\ -1.16 X 10"^") 



TABLE II: 1(7 limits on the parameter values for the beyond-detailed-balance scenario, for positive and negative curvature, and 
for two values of the effective neutrino species parameter AA^^ (see text). 








a 


A 


(-0.01,0.01) 


(0,2) 


(0, 1) 


(0.98, 1.01) 



1.1 



1.06 



1.02 



0.98 



0.94 



0.9, 





-6 -4 -2 2 4 
"ko x10'^ 

FIG. 4: (Color online) Contour plots of different pairs of 
free parameters in the hey and- detailed-balance scenario, un- 
der SNIa, BAO and CMB observational data. In each case 
the parameters not included in the plots have been marginal- 
ized over. Color scheme as in Fig. O 



As we observe, in Icr confidence the running parameter 
A of Hofava-Lifshitz gravity is restricted to the interval 
|A — 1| < 0.02, for the entire allowed range of ^3 (tliat is 
of 0-3). Finally, the best fit value for A restricts |A — 1| to 
much more smaller values, namely |Af,./ — 1| ss 0.002. 



TABLE III: 1(7 limits on the free parameters of the beyond- 
detailed-balance scenario. The cosmological parameters Q.mO, 
^bo, and Hq have been marginalized over. 



IV. PHASE-SPACE ANALYSIS OF 
HORAVA-LIFSHITZ COSMOLOGY 

In this section we review the results of the phase-space 
and stability analysis of Hofava-Lifshitz cosmology, with 
or without the detailed-balance condition, following . 
We are interested in investigating the possible late-time 
solutions, and in these solutions we calculate various ob- 
servable quantities, such are the dark-energy density and 
equation-of-state parameters. 

We start by transforming the cosmological equations 
into an autonomous dynamical system [45|, introducing 
suitable dimensionless variables which are combinations 
of the model variables and parameters. Then we extract 
the critical points of the autonomous system, and in order 
to determine their stability we linearize it around them 
and we examine the eigenvalues of the corresponding co- 
efficient matrix of the perturbation equations. 

In the case where the detailed-balance condition is im- 
posed, we find that the universe can reach a bouncing- 
oscillatory state at late times, in which dark-energy, be- 
having as a simple cosmological constant, will be domi- 
nant. Such solutions arise purely from the novel terms 
of Hofava-Lifshitz cosmology, and in particular the dark- 
radiation term proportional to is responsible for the 
bounce, while the cosmological constant term is respon- 
sible for the turnaround. 

In the case where the detailed-balance condition is 
abandoned, we find that the universe reaches an eternally 
expanding solution at late times, in which dark-energy, 
behaving like a cosmological constant, dominates com- 
pletely. Note that according to the initial conditions, the 
universe on its way to this late-time attractor can be an 
expanding one with non- negligible matter content, inde- 
pendently of the specific form of the dark-matter content. 
These features make this scenario a good candidate for 
the description of our universe, in consistency with ob- 
servations. Finally, in this case the universe has also a 
probability to reach an oscillatory solution at late times, 
if the initial conditions lie in its basin of attraction. 
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V. BOUNCE AND CYCLIC BEHAVIOR 

The possibility of late-time cyclic solutions that arose 
from the phase-space analysis, makes us to investigate it 
in more detail. Let us take a first look at how it is possible 
to obtain a cosmological bounce in this framework j46| . 
In the contracting phase we have H < 0, while in the 
expanding one we have H > 0, and by making use of the 
continuity equations it follows that at the bounce point 
H = 0. Throughout this transition H > 0. On the other 
hand, for the transition from expansion to contraction, 
that is for the cosmological turnaround, we have H > 
before and H < after, while exactly on the turnaround 
point we have H — 0. Throughout this transition H < 0. 

The above conditions for a bounce and a turnaround 
can be easily fulfilled in Hof ava-Lifshitz cosmology, as we 
observe from the two Friedmann equations (|10p and (1111) . 
In particular, a cyclic scenario could be straightforwardly 
obtained if we consider a negative dark radiation term 
and a negative cosmological constant. During the ex- 
pansion, the energy densities of all components decrease, 
which is not the case for the cosmological constant. Thus, 
its contribution will counterbalance that of dark matter, 
triggering a turnaround, after which the universe enters 
in the contracting phase. Then, after contraction to suf- 
ficiently small scale factors the dark radiation term will 
lead the universe to experience a bounce. Thus, the uni- 
verse in such a model indeed presents a cyclic behavior, 
with a bounce and a turnaround at each cycle [25j . 

The absence of singularities in a cosmological scenario 
is a significant advantage. However, one must examine 
the proceeding of fluctuations through the bounce. In 
general, non-relativistic gravities, such is Hofava-Lifshitz 
one, are usually able to recover Einstein's general relativ- 
ity as an emergent theory at low energy scales. Therefore, 
the cosmological fluctuations generated in this model 
should be consistent with those obtained in standard per- 
turbation theory in the IR limit [23]. In particular, the 
perturbation spectrum presents a scale-invariant profile, 
if the universe has undergone a matter-dominated con- 
tracting phase [13]. However, the non-relativistic correc- 
tions in the Hofava-Lifshitz action could lead to a modi- 
fication of the dispersion relations of perturbations. This 
issue has been addressed in [l^ (see [l^ and references 
therein for the perturbations of a pure expanding uni- 
verse in Hofava-Lifshitz cosmology), which shows that 
the spectrum in the UV regime may have a red tilt in 
a bouncing universe. Moreover, the perturbation modes 
cannot enter the UV regime in the scenario of matter- 
bounce. Thus, the analysis of the cosmological perturba- 
tions in the IR regime is quite reliable. 



VI. A MORE REALISTIC HORAVA-LIFSHITZ 
DARK ENERGY 



through (O,® in the detailed-balance case, or through 
([T ^ . dT^ in the beyond-detailed-balance case. Thus, one 
can straightforwardly obtain the dark-energy equation- 
of-state parameter in both cases, as wde = Pde/pde- 
As can be immediately seen, in both cases wde lies above 
the phantom divide. However, according to observations, 
Wde could have crossed —1 in the recent cosmological 
past. Therefore, the question is wether we can formulate 
an extension of Hofava-Lifshitz cosmology, in which the 
dark energy equation-of-state parameter can experience 
the phantom-divide crossing. 

For this shake we allow for an additional scalar field, 
which will contribute to the dark energy sector [s^] ^■ 
Hence, we add a second scalar ct, with action 



3A- 1 0-^ 
4 7V2 



-I- 



/ii/i2crVV - ^hlaV^cT + ^hjaV^a - V{(j) 



,(22) 



where V{a) accounts for the potential term of the cr-ficld 
and hi are constants. Assuming homogeneity, that is 
a = cr(i), its evolution equation will be given by 



ij + 3H& 



dV{a) 



3A - 1 da 



0. 



(23) 



Additionally, it can be easily seen that its contribution 
to the Friedmann equations of section |lT] will be the stan- 
dard scalar-field one, and thus one can absorb it in an 
extended dark energy sector, with energy density and 
pressure given by: 



Pde 



— - — (7 +V{a) + ^^^ TT-r + 



o 2 2 A 2 



Pde 



3A 



8(3A - l)a4 8(3A - 1) 

4 " 8(3A - l)a4 ~ 8(3A - 1) 

in the detailed-balance case, and by: 



Pde\, 



3A 



4 

3A - 1 



+ T/(a) + :^ + 



ai (J3K OiK 



V{a) 



+ 



6a4 



6a6 



in the beyond-detailed-balance one. Note that the dark 
energy density in both cases satisfies the usual conserva- 
tion equation. 

The aforementioned extended version of Hofava- 
Lifshitz dark energy can have a very interesting phe- 
nomenology. Firstly, the corresponding equation-of-state 
parameter wde can be above —1, below —1, or ex- 
perience the — 1-crossing during the cosmological evo- 
lution, as can be straightforwardly seen by the ratio 
Wde — Pde/ Pde- Thus, in this case, artifacts of Hofava- 
Lifshitz gravity could be detected through dark energy 



In section [H] we formulated Hofava-Lifshitz cosmology, 
in which one can define the effective dark energy sector 



^ Note that one could alternatively generalize the gravitational 
action of Hofava-Lifshitz gravity itself 0. |32|. |48|| . 
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observations. However, one still cannot distinguish be- 
tween this model and alternative models that allow for 
the realization of wde < — 1 phase, such are modified 
gravity or models with phantom jsoj or quintom 
fields j_51j. However, note that in the present formulation 
the additional scalar field is canonical, while in phantom 
and quintom scenarios the scalar field is phantom, and 
thus with ambiguous quantum behavior. The ability to 
describe the phantom phase and the phantom crossing 
with a canonical scalar field is a significant advantage of 
the scenario at hand, revealing the capabilities of Hofava- 
Lifshitz cosmology. 



VII. PERTURBATIVE INSTABILITIES IN 
HORAVA-LIFSHITZ GRAVITY 

In the previous sections we showed the advantages of 
Hof ava-Lifshitz cosmology at the background level. How- 
ever, despite the capabilities of the scenario, our analy- 
sis does not enlighten the discussion about the possible 
conceptual problems and instabilities of Hofava-Lifshitz 
gravity, nor it can address the questions concerning the 
validity of its theoretical background. Thus, in this sec- 
tion we are interested in performing a detailed inves- 
tigation of the gravitational perturbations of Hofava- 
Lifshitz gravity, using it as a tool to examine its con- 
sistency, studying both scalar and tensor sectors around 
a Minkowski background Q. 

We consider coordinate transformations of the form 
x'^ — ^ x'' = + Under this transformation the 
metric-perturbation around a given background changes 
as 6gfj,iy = 6gfj,iy - V^^^ - Vj,^^. Therefore, the general 
perturbations of the metric ([1]) read: 

dg,j = a^h.j - a^{d,Wj +djW,) -2a^^Pd^j +2a^d^djE. 

The vector modes are assumed to be transverse, that is 
diW^ = diQ^ = 0, while the tensor mode is forced to be 



transverse and traceless: dih^^ = 6^^ hi 



0. 



Let us now discuss the gauge fixing, which is required 
for the action derivation and the determination of the 
physical degrees of freedom. The projectability condition 
of Hofava gravity [3] requires that the perturbation of 
the lapse-function N depends only on time, thus (f) = 
(j){t). This allows us to "gauge away" the 4>- and B- 
perturbations, and also we can eliminate the Qi degree of 
freedom Q . Therefore, the remaining degrees of freedom 
are ^, E, Wi and hij. In summary, in the aforementioned 
gauge we obtain 

SN = 5N,^0 

= h^j - 2ijS,j + 2^^^JE - {^^WJ + ^JW^). (24) 

Note that since only perturbations imposed on the 
"same-time" spatial hypersurface are allowed, this is 
equivalent to a synchronous gauge choice. 



We now perturb the (prior to analytic continuation) 
Hofava-Lifshitz gravitational action up to second order. 
After non-trivial but straightforward calculations 0] , for 
the perturbed kinetic part of the action ^ we obtain 



^hijh'^ + (1 - 3A) (Stp^ - 2^V^E^ 



(1 - X)E\I^E 



(25) 



while for the perturbed potential part we acquire 

^2 



-—h-V^h'' 



8w^ 



2 2 a 

h ■w'^h'^ + ^ ^ h- W^h'^ 
32 + 32(1 -3A)'^ 



ajV(1-A)^,,^4 



4(1 -3A) 

27kVA^ 
16(1 - 3A) 

3^V^A^ 
'16(1 - 3A) 



4(1 - 3A) 



8(1 - 3A) 



EV^E 



(26) 



A. Scalar perturbations 

As can be observed from (pS)) , (|26p the action for scalar 
perturbations includes the two modes E and ip, and their 
equations of motion read: 

8" kV^(1-A)„2, kV^A , „ 
;^^+ 2(l-3A)^ -^+2(T^^ = ° (2^) 



8 1-3A 



-4'+— rrrz —^E 



4(1-A)(1-3A)^ 4(1-A)(1-3A) 



2(1 - 3A) 



' ^ ^ ^ = 0. (28) 



2(1 -3A) 



As can be seen these two equations are coupled, not al- 
lowing for a straightforward stability investigation. How- 
ever, we can still acquire information about the stability 
of the configuration by studying it at high and low mo- 
menta. Taking the IR limit of (|28|) . that is considering 
the low- A; behavior, it reduces to 



8 1 - 3A 



9«:V'A2 



k2 1- A ^ 4(1- A)(1-3A) 



= 0. (29) 



Thus, this decoupled equation acts as a low-momentum 
equation of motion for the scalar field ip. A straightfor- 
ward observation from (|29|) is that it leads to a ghost-like 
behavior, since it leads to the dispersion relation 



2 _ 2 

U! = m 



g^V^A^ 

32(1 - 3A)^ 



< 0, 



(30) 



which induces instabilities, regardless of the A-value and 
of the sign of the cosmological constant. Now, for high 
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fc, ([28]) reduces to 
8 1-3A 



V?/' = 0. 



k2 1 - A ' 2 (1 - 3A) 
Therefore, pil) yields the high-fc dispersion relation: 

1 - A 



9 n'^ip' 
J = — — 
16 



1 -3A 



(31) 



(32) 



B. Tensor perturbations 

Let us now examine the tensor perturbations. Their 
action can be extracted from (|25|) . ([26l) and therefore the 
graviton equation of motion writes as 



16(1 -3A) 



4w4 



16 



V^h'^ = 0. (33) 



Assuming graviton propagation along the direction, 
that is ki = = (0, 0, fc), the hij can be written as usual 
in terms of the Left and Right polarization components, 
and thus we derive the two equations for the different 
polarizations 



uj hL,_ 



^4, ,2 



— — fc hL.R± -r^k tiL.R 
16 iw"^ 



-^fc'Vfl = 0,(34) 



where the plus and minus branches correspond to Left- 
handed and Right-handed modes respectively. In this 
relation we have identified the light speed from the low k 
regime as = K'^fi^A/[16{l — 3A)]. The above equation 
system accepts a non-trivial solution only if the corre- 
sponding determinant is zero, which leads to the disper- 
sion relation 



c^k 



2 7„2 



4 2 4 4 

16 ^ Aw^ 



(35) 



C. Beyond Detailed Balance 



In order to avoid possible accidental artifacts of the 
detailed-balance condition, in this subsection we extend 
the investigation beyond detailed balance. As a demon- 
stration, and without loss of generality, we consider a 
detailed-balance-breaking term of the form Vii?jfeV*i?-''^. 
Thus, the corresponding contribution to the action will 
be 



5S, 



(2) 
bdb 



dt(fa 



6i/'VV 



(36) 



where rj is an additional parameter. It is straightforward 

(2) 

to calculate the modifications that SIJ^^ brings to the 



dispersion relations for scalar and tensor perturbations 
obtained above (expressions ([5^ and ([55)) respectively). 
The extended dispersion relations read: 



^2 (1 - A)' 



3«;2 (1 - A) 



16(1-3A)" 2(1-3A) 
for scalar perturbations (UV-behavior), and 



rjk^ 



16 



4^2 



4^ 



2 



(37) 



k"^ (38) 



for tensor perturbations. As was expected, the detailed- 
balanced-breaking term modifies mainly the UV regime 
of the theory. 



D. Instabilities 

Concerning the scalar perturbations, as was mentioned 
above ([29| . ([3T1) leads to instabilities. This unstable be- 
havior cannot be cured by simple tricks such as analytic 
continuation of the form /i — > i/x , —iw'^A [l3|, 

since in that case we straightforwardly see that the UV 
behavior is spoiled (see (|32|) ) and thus instabilities re- 
emerge at high energies. Even in this case though, we 
cannot evade the instability coming from the negative 
mass term, and thus IR instabilities persist as long as 
we have a non- vanishing cosmological constant. Finally, 
concerning the tensor sector, from psp we see that if we 
desire a well-behaved UV regime we cannot impose the 
analytic continuation. 

Proceeding to the relaxation of the detailed-balance 
condition, a crucial observation is that the ghost insta- 
bility of the scalar mode arises from the kinetic term 
of the action and thus the breaking of detailed balance, 
which affects the potential term, will not alter the afore- 
mentioned scalar-instabilities results. 



VIII. HEALTHY EXTENSIONS OF 
HORAVA-LIFSHITZ GRAVITY 

In the previous section we saw that Hofava-Lifshitz 
gravity in its simple version, with or without the detailed- 
balance version, suffers from instabilities and pathologies 
that cannot be cured. It is thus necessary to try to con- 
struct suitable extensions that are free of such problems. 

A quite general power-counting renormalizable action 
is M: 



with 



S — Skin + Si + S2 + Sn 



SMn^a J dtd''x^N[{K,jK'^~lK^)] 



(39) 



Si = / dtd-'x^N 



V9 
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Snew — J dt(fixy/gN [ai(aja') + 02(0^0*)^ + a^R^^aiaj 

+a4RS/ia' + a^\/iaj\/''-a^ + a6V*ai(aja-'') H ] . (40) 

Thus, apart from the known kinetic, detailed-balance 
and beyond-detailed-balance combinations that consti- 
tute the Hofava-Lifshitz gravitational action, in pOl) we 
have added a new combination, based on the term [9|]: 



N ' 



(41) 



which breaks the projectability condition, and the ellipsis 
in (|40|) refers to dimension six terms involving as well 
as curvatures. 

Such a new combination of terms seems to alleviate the 
problems of Horava-Lifshitz gravity, although there could 
still be some ambiguities [4l|. Therefore, one should re- 
peat all the investigations of the present work, for this 
extended version of the theory. 



IX. CONCLUSIONS 

In this work we reviewed some general aspects of 
Hofava-Lifshitz cosmology. Formulating the basic ver- 



sion of Hofava-Lifshitz gravity, with or without the 
detailed-balance condition, we extracted the cosmolog- 
ical equations. We used observational data in order to 
constrain the parameters of the theory, and amongst 
others we saw that the running parameter A, that de- 
termines the flow between the IR and the UV, is in- 
deed restricted in a very narrow window around its IR 
value 1. Through a phase-space analysis we extracted 
the late-time stable solutions, which are independent of 
the initial conditions, and we saw that eternal expan- 
sion, or bouncing and cyclic behavior, can arise naturally. 
Formulating the effective dark energy sector we showed 
that Hofava-Lifshitz cosmology can describe the phan- 
tom phase, without the use of a phantom field. However, 
performing a detailed perturbation analysis, we showed 
that Hofava-Lifshitz gravity in its basic version, suffers 
from instabilities. Thus, one should try to construct suit- 
able generalizations, that are free from pathologies, and 
then repeat all the above steps of cosmological analysis. 
Such a task proves to be hard, but it is necessary if we 
desire Hofava-Lifshitz gravity to be a candidate for the 
description of nature. 
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